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ABSTRACT 


The  ability  to  predict  torsional  instability  in  the 
early  stages  of  design  can  have  important  consequences 
on  the  design  of  both  conventional  and  high  performance 
ships.  This  thesis  develops  fast  approximate  methods  of 
torsional  buckling  analysis  for  particular  application 
in  the  concept,  feasibility,  and  preliminary  stages  of 
ship  design. 

Two  simplified  models  of  stiffeners  commonly  used  in 
ship  construction  were  presented.  The  first  was  an 
ideal  I  section  stiffener.  The  second,  a  flat  bar 
stiffener,  which  was  included  in  this  analysis  because 
of  oxits  favorable  productional  properties  and  wide 
uq^ge.  These  models,  with  initial  imperfections,  were 
subjected  to  axial  compressive  end  loading  and  the 
resultant  behavior  analyzed.  The  approach  of  the 
analysis  of  this  thesis  was  an  application  of  the 
energy  method  to  determine  the  critical  buckling  stress 
instead  of  the  more  commonly  used  equilibrium  approach. 
Both  beam  theory  and  thin  plate  theory  were  used  in 
energy-work  relationships  in  these  derivations.  Then 
the  first  yield  load  was  determined  utilizing  the 
Perry-Robertson  approach.  Integrated  into  all  phases  of 
this  exploration  was  the  concept  of  the  geometric 
imperfections  of  the  stiffeners. 
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NOMENCLATURE 


A  Nondimsnsional  paramster 

a  Length  of  stiffener  between  transverse  supports 
b  Uniform  stiffener  spacing 

b  Plating  effective  width 

e 

B  Angle  of  rotation  of  the  stiffener  about  the  toe 
* 

B  Initial  imperfection  of  the  stiffener  in  terms  of  a 
rotation  angle  about  the  toe 


B  Amplitude  of  the  sinusoidal  function  for  B 
o 

*  * 

B  Amplitude  of  the  sinusoidal  function  for  B 
o 

C  Rotational  spring  constant  (moment /length)  of  the 
supporting  plating 

C  Longitudinal  warping  constant  of  stiffener 
w 

D  Flexural  rigidity  of  plating 

0  Flexural  rigidity  of  web  plating  of  stiffener 
w 

d  Depth  of  flat  bar  stiffener;  overall  depth  of  Tee 
stiffener 

d  Depth  of  stiffener  to  midthickness  of  flange 
c 

d  Depth  of  web 
w 

E  Young's  Modulus  of  the  material 


H  Membrane  strain 

X 

f  Width  of  stiffener  flange 

M 

6  Shear  Modulus  of  the  material 

h  Height  of  neutral  axis  of  plate-stiffener 
combination  from  midplane  of  plating 

I  Effective  vertical  moment  of  inertia  of  stiffener 
and  associated  effective  width  of  plating 

I  Polar  moment  of  inertia  of  stiffener  about  toe 
P 

I  Vertical  moment  of  inertia  of  stiffener  alone  about 
t 

toe 

X  Moment  of  inertia  of  stiffener  about  web  plane 
z 

J  St.  Venant*s  torsion  const  an’;  for  stiffener 
L  Reduced  slenderness  ratio 

M(x)  Moment  developed  in  the  stiffener  by  the  total 
torque 

m  Mode  number 
P  Axial  end  load 

s  Haight  of  stiffener  shear  center  above  toe 

ft  Elastic  axial  torsional  buckling  stress  in 
CR 

stiffener  (critical  stress) 


ft  Elastic  axial  torsional  buckling  stress  limit  of 
CL 

stiffener  (classical  stress) 


ft  Axial  stress  in  stiffener 

e 
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ft  Axial  stress  component 


X 


(ft  ) 

X  MAXT 


Maximum  compressive  stress  developed 
the  torque 


by 


ft  Tensile  yield  stress  of  material 

y 


T  Amplitude  of  the  sinsoidal  torque  function 
o 


T  Total  torque  developed  in  the  stiffener  structure 
by  the  loading 

T  St.  Venant*s  torque  developed  in  the  stiffener 
SV 


T  Uarping  torque  developed  in  the  stiffener 

M 


t  Plate  thickness  m 

t  Stiffener  flange  thickness 
f 

t  Stiffener  web  thickness 
w 

u  Displacement  in  the  longitudinal  (x)  direction 

V  Total  strain  energy  of  the  structure 

V  Sideways  flexure  (y)  of  stiffener 
* 

V  Initial  imperfection  of  stiffener  in  horizontal  (y) 

direction 

U  Total  work  of  the  applied  force 

w  Vertical  flexure  (z>  of  stiffener 
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« 

w  Initial  imparf act  ion  of  the  stiffener  in  the 
vertical  (z)  direction 

z  Height  of  stiffener  centroid  above  toe 

^  Poisson’s  Ratio 


Note:  K,y  or  z  subscripts  indicate  partial  derivatives 

with  respect  to  those  coordinates  .  These  are  used 

*  *  * 
particularly  with  v,  v  ,  w,  w  «  B,  B  . 
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INTRODUCTION 


The  subject  of  this  thesis  is  to  develop  a 
simplified  analytical  method  of  predicting  torsional 
buckling  of  stiffeners  used  in  ship  structures. 
Torsional  buckling  or  instability  is  the  phenomena  by 
which  a  column,  at  a  certain  axial  load,  fails  suddenly 
in  a  combined  mode  of  twist  and  lateral  bending  of  the 
cross  section  (reference  7>.  Figure  (0-la)  (from 
reference  6)  illustrates  the  failure  mode  referred  to 
here  as  torsional  buckling  -  in  this  case  for  a  flat 
bar  stiffener.  Figure  (O-lb)  illustrates  the  failure 
mode  of  the  tripping  phenomenon  which  is  not  covered  in 
this  analysis.  Figure  (0-2)  (reference  6)  is  the 
experimental  load  deflection  curve  for  the  failure  of  a 
flat  bar  stiffener  under  an  axial  load.  Point  P  of  the 
solid  (experimental)  curve  of  figure  (0-2)  is  the 
ultimate  strength  load.  As  can  be  seen  from  figure  (0- 
2),  the  curve  proceeds  in  a  downward  fashion  -  the  bar 
looses  stiffness  and  fails  from  this  point. 

The  ability  to  predict  torsional  instability 
in  the  early  stages  of  design  can  have  important 
consequences  on  the  design  of  both  conventional  and 
high  performance  ships.  It  is  the  intention  of  this 
thesis  to  develop  fast  approximate  methods  of  torsional 
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Figure  (0-la)  Torsional  Buckling 


Figure  (0-lb)  Tripping 


Figure  (O-l)  Characterization  of  Torsional  Buckling 
and  Tripping  in  a  Flat  Bar  Stiffener 
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If.  H  a  ut 


£  oa 


Figur* 


(0-2)  Experimental  Load  Deflection  Curve 
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buckling  analysis  fon  particular  application  in  the 
concept,  feasibility,  and  preliminary  stages  of  ship 
design. 

To  this  end,  two  simplified  models  of  stiffeners 
commonly  used  in  ship  construction  will  be  presented. 
The  first  is  an  ideal  I  section  stiffener.  The  second, 
a  flat  bar  stiffener,  which  is  included  in  this 
analysis  because  of  its  favorable  product ional 
properties  and  wide  useage.  These  models,  with  initial 
imperfections,  will  be  subjected  to  anial  compressive 
end  loading  and  the  resultant  behavior  analyzed.  The 
approach  of  the  analysis  of  this  thesis  will  be  an 
application  of  energy  methods  to  determine  the  critical 
buckling  stress  instead  of  the  more  commonly  used 
equilibrium  approach.  Both  beam  theory  and  thin  plate 
theory  are  used  in  energy-work  relationships  in  these 
derivations.  Then  the  first  yield  load  is  determined 
utilizing  the  Perry-Robertson  approach  as  described  in 
references  (9,14).  Integrated  into  all  phases  of  this 
exploration  will  be  the  concept  of  the  geometric 
imperfections  of  the  stiffeners. 

It  is  the  intention  of  this  thesis,  using  the  above 
mentioned  approach,  to  lay  groundwork  which  would  have 
the  potential  for  future  expansion  to  include  mode 
interaction  -  to  be  able  to  easily  analyze  combined 
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loadings  etc. . 

Particular  emphasis  is  placed  here  upon  inclusion 
of  the  geometrical  initial  deformities  of  the 
stiffeners  Mhich  occur  during  the  manufacture  and 
through  day  to  day  usage  of  the  ship.  It  is  well  known 
that  these  geometric  imperfections  can  drastically 
reduce  the  strength  of  the  structural  member.  But 
current  design  standards  (references  10,11  etc)  have 
not  tried  to  quantify  these  values.  Some  design  codes 
(example:  reference  10)  have  specified  tolerence  limits 
of  various  geometric  parameters  which  must  be  met 
before  the  design  ,  formulae  are  considered  valid. 
Granted,  in  the  concept,  feasibility  and  preliminary 
stages  of  ship  design  where  the  ship  only  exists  on 
paper,  the  designer  would  have  no  knowledge  of  the  end 
resulting  imperfections  of  the  actual  ship.  Hopefully, 
in  the  future,  with  the  advent  of  better  production 
control  in  shipyards  and  the  incursion  of  the  computer 
into  ship  production  facilities,  this  information  will 
be  collected  and  statistically  analyzed. 

In  addition,  once  the  ship  is  manufactured,  these 
formulae  will  provide  a  quick  and  easy  evaluation  of 
the  actual  critical  buckling  stress  in  terms  of 
physical  parameters.  Thus,  this  analysis  has  made  a 
special  effort  to  define  the  imperfections  in  easily 
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maasured  terms.  This  analysis  measures  the  horizontal 
or  vertical  movement  of  the  shear  center  of  the 
stiffener  at  its  point  of  greatest  deflection  and 
assumes  a  simple  sinusoidal  distribution.  Knowledge  of 
the  web  depth  easily  converts  the  imperfection  into  a 
rotational  angle.  Here,  the  stiffener  is  considered  to 
rotate  about  its  toe  or  base.  In  other  words,  the  base 
plating  is  considered  much  larger  and  stiffen.  Figure 
(0-3)  illustrates  this  simplified  concept  for  a  Tee 
stiffener. 


Figure  (0-3)  Tee  Stiffener 

Currently  the  method  of  handling  the  stiffener 
failure  problem  is  incorporated  into  the  use  of  factors 
of  safety.  But,  with  recent  attempts  to  further 
optimize  ship  structures  -  in  particular  to  minimize 
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tha  amount  of  weight  of  a  ship  that  must  be  devoted  to 
its  structure  -  better,  yet  reliable,  methods  must  be 


devised.  Remember,  with  a  ship,  less  weight  in  the 
structure  equates  to  more  weight  for  the  payload  be  it 
weapons  for  a  military  ship  or  cargo  for  a  civilian 
ship. 

This  thesis  concludes  with  a  comparison  between  the 
formulae  derived  from  the  simplified  models  presented 
here  and  published  formulae,  published  finite  element 
analysis,  and  experimentally  derived  results. 
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Figura  (1-lb)  Undeformad 


Figura  (1-lc)  Daformad 


Figura  (1-1)  Charactarization  of  Torsional  Buckling 
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1,1  Model  la. 

The  first  model  considered  is  a  simplification  of 
an  ideal  I—section.  It  consists  of  an  I-beam  with  zero 
torsional  rigidity.  The  web  of  the  I-beam  transmits 
only  shear  forces.  In  effect,  it  looks  like  two  unequal 
flanges: 


I-beam  Stress  Distribution 

Figure  <l-2)  Model  I 

This  model  has  been  considered  in  reference  (4). 
However,  reference  (4)  used  an  equilibrium  approach  and 
initial  imperfections  were  considered  as  an 
eccentricity  of  the  axial  load.  The  resulting  equations 
of  the  critical  stress  were  quite  complicated. 

Figure  (1-3)  illustrates  the  geometric  parameters 
of  this  model. 
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Figurs  (1-3)  -  Geometric  Torsional  Buckling 

Parameters  for  Model  I 


-  Moment  of  inertia  about  the  web  plane  (top  flange 

only) 

3 

I  ■  (t  f  )/12 

2  f  W 

-  Height  of  shear  center  above  too  (origin)  (top 

flange  only) 


c 

-  Longitudinal  warping  constant  about  shear  center 
of  stiffener  alone  (top  flange  only) 

C  *0 


-  Vertical  moment  of  inertia  about  toe  (top  flange 
only) 

2  2 
X  -  f  t  (d  •••  (l/12)t  ) 


I  -  Polar  momant  of  inertia  about  toe 
P 

I  -  I  +  I 
p  t  z 

z  -  Height  of  centroid  above  toe 
z  «  d 

c 

J  ~  St.  Venant*s  torsion  constant  (top  flange  only) 
3 

J  ■  (ft  )/3 
w  f 
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1.2  Modal  II. 

Tha  second  modal  considered  is  the  simple  flat  bar 
stiffener  illustrated  in  figure  (1-4). 


Figure  (1-4)  Model  II  -  Flat  bar  stiffener. 

Figure  (1-5)  illustrates  the  geometric  parameters 


of  this  model 


i. 


I 

B 


i 

i 


» 


L. 


Figure  (1-5)  -  Geometric  Torsional  Buckling 

Parameters  for  Model  II 


I  -  Moment  of  inertia  about  the  meb  plane 
2 

(web  only) 

3 

I  »  <d  t  )/12 

2  W  W 

s  •  Height  of  shear  center  above  J^oe  (origin) 

(web  only) 

B  »  d  /2 

w 

C  Longitudinal  warping  constant  about  shear 

w 

center  of  stiffener  alone 

C  -  0 
w 

I  -  Vertical  moment  of  inertia  about  toe  (web  only) 
t  3 

I  «  (t  d  )/3 
t  w  w 

I  *>  Polar  moment  of  inertia  about  toe 
P 

I  ■  I  ♦  I 
P  t  2 
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z  ~  Height  of  centroid  above  toe 


z  »  d  /a 

w 

J  -  St.  Venant*s  torsion  constant 
3 

J  -  <d  t  )/3 
w  w 

D  -  Flexural  rigidity  of  the  web  plate 
w 

3 

Et 

w 

w  a 

13(1-  V  ) 
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2.  MODEL  I  -  TORSIONAL  BUCKLING  UNDER  AXIAL  LOADS 


2. 1  Development  of  tbS  Equation. 

Since  Model  I  basically  consists  of  two  unequal 
flanges  with  the  top  flange  notating  about  the  bottom 
flange,  the  torsional  rigidity  is  considered  only  for 
the  top  flange.  Thus  the  strain  energy  equation 
consists  of  three  parts  (reference  1).  The  first  term 
represents  sideways  bending.  The  second  term  represents 
longitudinal  warping.  The  third  term  represents  the 
torsional  rigidity  of  the  top  flange  alone.  The  strain 
energy  for  a  length  a  is: 

^a  2  2  2 

V  »  1/2  J  <EI  V  +  EC  B  +  GJB  )dK 

''0  Z  MX  W  XX  X 


Also  from  geometrical  properties  we  have  (see  figure 


1-3) : 

V  »  sB 

Thus: 

V  »  1/2 

V  *  1/2 


2  2  2  2 
(El  s  B  EC  B  -I-  GJB  )  dx 

Z  XX  W  XX  X 


2  2  2 
Ed  s  -»•  C  )B  GJB  I  dx 

Z  W  XX  X 


2S 


to 


B 

I 

I 


The  following  choice  of  B  is  made, 
meet  certain  boundary  conditions  (ie  B 
x^a  )  and  to  keep  B  simple* 


XH 


in  order 
30  at  x^O, 


B  “  B  sin  (raT^'x/a) 
o 

Thus  for  Model  I: 

2  2  2 
V  3  1/4  a  B  (m^'/a)  C  Ed  s 
o  z 


2 

+  C  )  (m  Tr/ak) 
w 


+  SJ  3 


i 

i 


It 
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2«  2  SiSD  Convent ipn 


Before  development  of  the  virtual  work  equation,  a 
note  should  be  said  about  the  sign  convention  used 
here.  Compression  is  held  to  be  a  positive  stress  and 
tension  as  negative.  The  following  figures  illustrate 
the  sign  convention  used  here. 


Figure  (2-1)  Sign  convention  for  inplane  axial  and 

shear  stress  components. 
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2, 3  DgvgioglSSDt  of  yiCtiifiii  ySllE  ESiiStiSDa. 

Modsl  I  is  not  assumed  to  be  perfect  but  rather  to 

have  initial  imperfections.  These  initial  imperfect  ions 

*  «  « 

V  ,M  and  B  and  the  additional  deformations  v,w  and  B 
are  illustrated  in  figure  (2“3). 


Figure  (2-3)  Illustration  of  initial  imperfect  ions 

and  additional  deformations  B. 


The  virtual  work  for  the  case  of  an  inplane  anial 
load  becomes  (reference  l)i 


JJ"  %  u(y, z)  dy  dz 


A  M 

Thus  using  figure  (1-1),  it  can  be  shown  that  the 
strain  of  the  centroidal  an  is  in  the  bent 
configuration,  •  (membrane  strain),  is  (reference  l)i 

M 


2  Z 

•  ■  -u  l/2(v  )  *  i/2(w  ) 

K  K  H  X 
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for  a  parfect  I  saction.  With  im par fact  ions, 

subst itutai 
* 

V  +  V  for  V 

« 

M  -tt  w  for  w 


Thus* 

*  2  *  S 

*  ■  -u  l/2<v  +v  )  +  l/2(w  ) 

XX  K  X  K  X 

using  tha  inaxtansibility  assumptions,  sat: 

Ir  «  0 

X 

wa  hava: 

«  2  *  2 
u  ■  l/2<v  +v  )  +  l/2(w  +w  ) 

X  XX  XX 


r: 


■a  *  2  *2 

u  ■  1/2  I  C(v  +v  )  +  (w  ■♦•w  >  3  dx 

XX  XX 


Mhara  u  »  u<y, z). 


From  geomatrical  ralat ionships  (saa  figura  (2-3)) 


V 

a 

zB 

w 

ffi 

>* 

1 

a 

» 

* 

» 

* 

V 

a 

zB 

w 

a  -yB 

IX) 

1 

V 

a 

zB 

w 

»  -yB 

X 

X 

X 

X 

« 

* 

« 

* 

V 

a 

zB 

w 

a  -yB 

X 

X 

X 

X 

Substituting  and  intagrating  wa  get 


N  mi  ij  ■.  j  ^ 


/-«  2  2  2  *  2  *2 
u(y,  z)  •  1/2  f  (zB  ■t>2zBB  -t-z  B 

^0  M  X  K  X 


2  2  2  *2  2  *2 
+yB  +2yBB  +yB  )dx 

X  XX  X 


✓a  2  e  *  c 

u(y,z)  «  1/2  j  (z  +  y  ) (B  +  B  )  dx 

XX 


2  *  2 
•  (B 

X  X 


Again  using  the  simple  choice  for  B  and  B  : 


B  *  B  sin(ffl  t'  x/a) 
o 

*  * 

B  «  B  sin(m'^x/a) 
o 


Thuss 


u<y,z) 


■  1/2 


2  2 

(z  ■f'y  )C(m/^/a)B  cos(m7)’x/a> 

o 


*  2 

*  ((m^/a)B  cos(m  r^x/a)  3  dx 
o 


u<y,z)  «  1/2 


2  2  *  2 
(z  *y  )  (B  -t-B  ) 
o  o 


2 

*  C<m 7r/a)cos(ffl ^x/a) 3 


dx 


2  2  *2  2 
u(y,z)  •  l/2(z  *y  )  (B  i-B  )  (mT^/a) 

o  o 


*  Cx/2  ♦  (sin(2m7^x/a)/(4m;^/a>3 

0 
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u(y, 2) 
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2.4  Dsiscmicatign  s£  ibS  QClilSSi  SUSBlica  iiCSSSs. 

Applying  th«  calculus  of  variations  with  respsct 

to  B  to  the  work  and  strain  energy  equations  yields 
o 

(reference  S) i 

*  2  2 

<fu  •  1/2  S  <B  +B  >  {m?^'/a)  a  I  Xb 
e  o  o  P  o 

2  2  2 
(Tv  -  1/2  a  B  (niD'/a)  CE(I  s  +  C  >  (m /V/a) 


6J3  a  B 

o 

Applying  the  principle  of  minimum  potential  energy* 

TT  m  Sv  ~  Su  •  0 


we  get! 

(5^  M  * 

Solving  for  *  : 


2  _  2 
CE(I  s  C  )(m///a)  +  GJ3 

z  w 


B  B 


(2-2) 


It  can  be  seen  by  inspection  that  the  lowest 
buckling  stress  occurs  for  one  wave,  m  »  I,  since  m 
must  be  an  integer.  Thus* 
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£m  S  Determinat  ion  s£  &t!S  EiCSl  Yi.gld  (eSSMs. 

To  datemiinB  th®  fii^t  yield  load*  the  total  torque 
developed  in  the  torsional  buckling  process  must  be 
first  evaluated.  This  torque  will  consist  of  two  parts. 
The  first  term  consists  of  the  well  known  St.  Venant’s 
Torque  (reference  1).  The  second  term  represents  the 
warping  torque  (reference  2). 

From  reference  ( 1 ) s 

T  »  T  +  T 

SV  W 

where! 

T  «  GJB 
SV  K 

and  from  reforonce  (2) i 
T  -  -  El  d  w 

U  Z  C  MXX 

and  from  geometry; 
w  -  d  B 

M  C  K 

resulting  ini 

2 

T  ■  -  El  d  B 

U  z  e  HHX 

Thus  the  total  torque  isi 
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(2-5) 


2 

T  «  6JB  -  El  d  B 

K  Z  C  MXX 

Figura  (2-3)  illustrates  the  above  equations. 


Again  a  simple  choice  of  B  is  made,  and  with  m  «  1: 

B  »  B  sin(7^x/a)  (2-6) 

o 

Equation  (2-5)  now  becomes! 

T  ■  SJ(i^/a)B  cosOT-x/a) 
o 

2  3 

El  d  (  TT’/a)  B  co8(^^x/a) 

2  c  o 


2  2 

T  *  ()f-/a)B  cos(7^x/a)CSJ  +  El  d  (7^/a)  3 


o  z  c 

Choosing  a  simple  sinusoidal  shape  for  T  allows: 

T  -  T  cos  (7?  x/a) 
o 

Thus! 


T 


(77’/a)B  CBJ  +  El  d  (77' /a) 


2 

3 


(2-7) 


o 


o  z  c 


To  eliminate  B  ,  equation  (2-4)  is  applied  to  equation 
o 

(2-7). 

* 

B 

o 

2  2 

T  ■  CBJ  +  El  d  (7J^/a)  3(^/a) - 

o  z  c 

(«  /«  -  1) 

CL  CR 

(2-8) 
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From  »Bctional  considerations,  an  evaluation  of  the 
maKimum  stress  can  be  made. 


5 


\ 


Figure ^<2-5)  Sectional  illustration. 

The  moment  developed  by  the  torque  in  terms  of 
deflections  is  <reference  2) i 
M<k)  ■  -  EI  n 

2  KM 

So  the  maKimum  compressive  stress  developed  in  the 
flange  due  to  the  torque  ist 


-  EI  d  B  f 

Z  C  KK  W 


(k  } 

M  MOXT 


21 


(ft  }  -  -  Ed  B  f  /2 

M  MflXT  C  KK  w 


Application  of  the  above  definition  for  B  (equation 
2-4)  leavesi 
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I. 


Ed  f 

c  w  2 

(S  )  - - i.T^/ai'i  B  sin<77'x/a) 

X  MflXT  o 

2 

A  choica  of  x  is  made  such  that  <9r  )  is  a  maximum: 

X  MAXT 


Ed  f 

c  w  2 

(K  )  - - (  TT'/a)  B  (2-9) 

X  MAXT  o 

2 


Equations  (2-7)  and  (2-9)  are  combined  to  eliminate  B 


Ed  f  T  (  /a) 

c  w  o 

(S  )  - - 

X  MAXT  2  2 

2CSJ  +E1  d  (7»^/a)  : 
z  c 


(2-10) 


Using  equation  (2-8)  leads  to: 


Ed  f  (JT'/a) 

C  M 


(ft  ) 

X  MAXT 


B 


Eft  /ft  -  13 
CL  CR 


(2-11) 


Note  that  when  B  «  0  that  this  equation  leads  to  an 

o 

indeterminate  form  of  ft  since  ft  «  ft  . 

X  CL  CR 

The  total  maximum  stress  can  now  be  described  by: 


(ft  )  «  (ft  )  +  P/A 

X  TOTAL  X  MAXT 


where 
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P/A  «  « 

CR 


Thus  ths  yield  stress  has  been  shown  to  be: 

ft  «  (ft  )  ft  (2-12) 

y  M  MAXT  CR 

Using  equation  (2-11)  in  equation  (2-12): 

2  * 

Ed  f  (7?" /a)  B 


c  w 


2  ft  /ft  -  1 

CL  CR 


2 

Solving  for  ft  and  dividing  by  ft  to 
CR  y 

nondimensional  form  yields: 


(2-13) 

place  into 


(2-14) 


Choosing  ths  smallest  root  and  simplifying  the  equation 
by  introducing  a  slenderness  ratio  L: 


(2-15) 
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and  lattingj 


2  * 

<7r/a)  Ed  f  B 

C  M  O 


(2-lS) 


Then  : 


(1  +  1/L  +  fl)  ~  \  <1  +  i/L  +  A)  -  4/L 


<2-17) 

A  plot  of  ft  /ft  vs  L  for  various  values  of  A  is 
CR  y 

shown  in  figure  (2-6).  Figure  (2-6)  was  compared  with 
various  design  criteria  from  reference  (10)  and  it  was 
found  that  equation  (2-17)  was  similar  to  design 
equations  of  reference  (10).  However,  the  reduced 
slenderness  ratio,  L,  of  reference  (10)  is  based  upon 
simple  Euler  buckling  and  the  geometric  parameter.  A, 
is  of  course  different.  The  column  curves  of  reference 


(10)  were  quite  similar  in  appearance  and  range  of 
values  to  that  found  in  figure  (2-6). 
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3.  MODEL  II  -  TORSIONAL  BUCKLING  UNDER  AXIAL  LOADS 


I  Deygloowent  of  the  Str§i.n  Enerfl^  Equation. 

Since  Model  II  is  the  very  simple  case  of  a  flat 
bar  stiffener,  thin  plate  theory  can  be  easily  utilized 
in  the  development  of  the  strain  energy  equation.  Much 
of  this  section  has  been  patterned  after  a  similar 
analysis  presented  in  reference  (3). 

For  a  web  plate  of  uniform  thickness  and  a  length  a 
the  strain  energy  is  repensented  by  (reference  5) : 


making  a  simple  choice  for  B; 

B  ■  B  sin  (m^x/a) 
o 

substituting  into  the  expression  for  the  strain  energy 
and  integrating  over  a  constant  depth  d  results  in: 
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3.  £  Dsvs^02msQt  of  %hB  V irt yal,  Work  Equat  ion,, 

As  for  Modal  1,  Modal  II  is  not  assumed  to  be 

perfect  but  rather  to  have  initial  imperfections.  These 

*  «  * 

initial  imperfections  v  ,  w  ,  6  and  the  additional 

deformations  v,  m,  and  B  are  illustrated  in  figure 
(3-1). 


Figure  (3**1)  Illustration  of  the  initial  imperfect  ions 
V*, m*,B*  and  additional  deformations  v, w, B. 


For  the  case  of  an  inplane  axial  load,  the  work 
equation  is  i 


W 


%  u<ytz) 
K 


dy  dz 


From  figure  (3-1) «  it  can  be  shown  that  i 

£  £ 

-  u  ♦  l/£(v  )  l/£<w  ) 

MM  X 


•  - 

M 
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for  a  perfect  web  section.  With  imperfections,  and 
using  the  inextensibility  assumption  with  regard  to  the 
length  of  the  stiffener  setting  B  *  ®»  the  result  is  : 


«  2  *  a 

u  =  1/2  (v  +  V  )  +  l/2(w  +  w  ) 

M  XX  XX 


u  *  1/2 


»  2  *2 
C(v  V  )  +  l/2(w  w  >  3  dx 


o  X  X 


X  X 


(3-2) 


where  u  «  u(y, z) 

As  shown  in  figure  (3-1),  the  geometrical 
relationships  (equations  2-1)  are  used  to  simplify 


equation  (3-2) i 


u(yyz)  *  1/2 


'a  2  2  * 

(z  +  y  )  (B  B  )  dx 


A  simple  choice  is  made  for  B  and  B 
B  ■  B  sin  (mT^x/a) 


*  » 

B  »  B  sin  (m/T^x/a) 


Thus: 

2  2  *  2 

u<y, z)  ■  l/4(z  y  )  (B  +  B  >(mJ!^/a)  a 

o  o 

The  work  equation  now  becomes: 

rr  2  2*2 

W  >  J  /  *  Cl/4(z  +  y  )  (B  B  )  (m:?'/a>  a3  dy  dz 
^'-'A  X  o  o 


o  o 
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*  2  2 

W  >  1/4  K  (B  -«■  B  )  (M^/a)  a  I 

X  o  o  I 


3.3  Diytemiiinat ion  gf  tha  Critical  BuchliDfl  Stress. 

Applying  the  calculus  of  variations  with  respect  to 

B  to  the  Mork  and  strain  energy  equations  yields: 
o 

-  *  2  2  r 

dU  >  1/2  B  <B  B  >  (m/'^/a)  al  d  B 


K  O  O 


P  o 


.  2  2  2 
a  V  ■  <  D  ad  <M^/a)  B  C  d  (m  7^/a) 

M  O 

+  6(1  -V)  3/6  +  CB  a/2  >Xb 


Applying  the  principle  of  miniMum  potential  energy: 

7T  *  i'v  -  Sw  «  0 
Su  • 

Solving  for  Baft  : 

e  CR 

2  2  2 

D  d(m7^/a)  Cd  (m  7^/a)  +  6(l-V)3/3  +  C  B 


(m  /a) 


B  B 
o  o 

(3-3) 


For  the  case  of  zero  restraint  against  rotation) 


C  »  0: 
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2  2  2 

D  d<in7r/a>  Cd  (in>/a)  +6(1-V)3 

M 

ft  - - 

CR 

3  I 

P 


B 

o 

* 

B  +  B 
o  o 

(3-4) 


It  is  sasily  sesn  that  the  lowest  ft  for  this  case 

CR 

results  from  the  simplest  of  mode  shapest  m=l. 

When  the  rotational  restraint  is  not  equal  to  zero; 

the  lowest  value  of  ft  depends  upon  the  degree  of 

CR 

rotational  restraint: 


3  1/4 

m  *  (a/7^  )  <3C/D  d  ) 

w 


(3-5) 


With  the  use  of  equation  (3-5),  m  may  be  evaluated 
and  rounded  up  or  down  to  the  nearest  integer  for  the 
correct  wave  shape. 

It  is  easily  seen  that  as  B  becomes  large  ft 


CR 


approaches  a  limit.  Thus  with  (ft  ) 

CR  L.TMIT 


ft  I 

CL 


ft 

CL 


2  2  2 

0  dCmT^Va)  Cd  (m7^/a>  +  6(l-V)3/3  +  C 

w 


2 

(m  "J^/a)  I 

P 


(3-6) 


Figure  (3-2)  illustrates  the  behavior  of  ft  with 

#  CR 

increasing  B  B  . 

Further  simplification  and  rearranging  of 


•quations  (S-S)  and  (3-6)  providesi 

« 

B 

o  ^  ^ 

o 

«  /ts  -  1 
CL  CR 


Figure  (3-2)  Behavior  of  K  with  B  *  B 

CR  o  o 


The  first  yield  load  Mill  not  be  determined  for 
this  case  under  this  method  of  analysis  due  to  the 
relative  complexity  of  equation  (3-3)  and  due  to  the 
fact  that  the  Perry-Robertson  method  Mill  not  lend 
itself  to  further  simplification  along  this  train  of 
analysis.  Since  a  simplified  formula  for  design 
purposes  is  the  goal  of  this  thesis,  the  next  chapter 
Mill  use  a  different  tactic  to  solve  the  problem. 
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4.  MODEL  II  -  6ENER;)LIZED  ANALYSIS 


4. 1  M»th9d  of  Approach  ^nd  §train  Ensrs^  Equat iqn. 

Tha  analysis  thus  far  for  ths  cass  of  a  flat  bar 
stiffensr  undsr  and  loading  has  procaeded  along  the 
linas  of  thin  plata  thaory  up  to  the  development  of  an 
expression  for  the  critical  buckling  stress. 

This  chapter  uses  beam  theory  to  develop  a  simpler 
formulation  in  a  manner  analagous  to  that  performed  on 
Model  1  in  chapter  2.  In  fact,  this  analysis  will 
merely  be  a  generalization  of  the  results  for  Model  I 
applied  to  Model  II.  It  must  be  remembered  here  that 
Modal  I  consistad  of  two  unequal  flanges  with  the  top 
flange  deforming  or  rotating  about  the  bottom  flange. 
At  no  time  in  the  derivation  of  the  formula  for  the 
critical  stress  of  the  tea  stiffener  is  there  an 
assumption  made  about  the  position  of  the  top  flange 
relative  to  the  bottom  flange.  Thus  the  same  analysis 
with  minor  changes  will  be  performed  for  the  flat  bar 
stiffener.  However,  in  the  case  of  the  flat  bar 
stiffener  of  Model  II,  there  is  a  restraining  effect 
provided  by  the  juncture  of  the  flat  bar  with  the 
plate.  This  in  affect  can  be  modeled  as  a  spring  with  a 
rotational  constant  C.  Thus  the  strain  energy  equation 
would  have  four  terms  (instead  of  the  three  terms  for 


Model  I>.  These  terms  represent  (in  order  of 
appearance)  sideways  bending;  longitudinal  warping; 
torsion;  and  rotation  of  the  supporting  plate  structure 
modeled  as  an  elastic  spring.  Thus  the  strain  energy 
equation  is: 


1/a 


f 


2  2  2  2 
CEI  V  •t'  EC  B  *  6JB  -f  CB  3  dx 

O  Z  XX  W  XX  X 


Again  choose: 

B  B  sin(m  H'x/a) 
o 

Substituting  and  integrating  results  in: 

2  2  4  2 

V  *  1/4  a  B  C  Ed  s  +  C  XmTT/a)  +  SJCrnTT/a)  +  C  3 
o  z  w 

Applying  the  calculus  of  variations  with  respect  to 

B  yields: 
o 

r  2  4 

a  V  *  1/2  a  B  C  Ed  s  +  C  )  (m  3^/a) 


GJ(m 


7^/a>^  +  C  3  X 


B 


(4-1) 


^•8  Ihs  iansilsQa. 

The  assumptions  used  in  deriving  the  work  equation 
for  Model  Z  are  the  same  as  for  Model  II.  Thus  from 
chapter  2  we  have: 

*  2  2 

M  *  1/4  «  (B  +  B  )  <m/^/a)  a  I 
e  o  o  p 


SI 
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reveals  the  only  difference  is  that  equation  (4-3)  has 
an  additional  term  in  the  numerator  (specifically  the 
spring  constant  C) . 

Since  it  is  very  difficult  to  quantify  the  spring 
constant  C  and  (hopefully)  since  its  value  is  small  in 
comparison  to  other  terms  in  expression  (4-3),  C  will 
be  set  to  zero  value  and  ignored  for  the  remainder  of 
this  analysis. 

Thus  for  the  case  of  C  «  0,  equation  (4-3)  reduces 
to  equation  (2-S). 

4.4  Determination  of  th§  F^rst  YiSi^  Load. 

The  determination  of  first  yield  load  for  the 
generalized  case  proceeds  in  a  similar  fashion  as  for 
the  case  of  Model  I. 

Again  the  total  torque  consists  of  two  parts 
(reference  1),  St.  Venant’s  torque  and  a  warping 
torque.  However,  in  this  case  the  warping  torque  is 
(reference  S) i 

T  ■  -El  d  w 

W  Z  KXX 

and 

w  ■  dB 

X  X 

resulting  ins 
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(4-8) 


Using  squat  ion  (4-8)  to  eliminats  B  i 

o 

2 

Ed  T  (T^/a) 
o 

K  MAXT  2  2 

CQJ  +  El  d  (^/a)  3 
z 


Using  aquation  (4-7)  in  aquation  (4-8)  laads  to: 

« 

B 

2  2  o 

(k  )  -  Ed  (TT'/a)  -  (4-9) 

X  MflXT 

(k  /k  -  1) 

CL  CR 

Again,  from  sactional  considerations,  the  total 
maximum  stress  can  be  described  by: 


k  >  (k  )  >  k 

y  X  MAXT  CR 


(4-10) 


Using  aquation  (4-9)  in  equation  (4-10),  solving  for 

2 

k  and  dividing  by  k  to  place  into  nondimensional 
CR  y 

form  yields: 

2  2  2  * 
k  Ed  (7>/a)  B 


i 


2  2  * 
(  /^/a)  Ed  B 


% 

y 


(4-12) 


Results  in  (equation  2-17) s 


2 

3  (1  +  1/L  +  fl) 

CR 


Ni 


2  2  2 
(1  +  1/L  fl)  -  4/L 


ft  2 

y 


(4-13) 


5.  COMPARISONS  UITH  OTHER  STUDIES 


5. 1  QsiSfiiCiSSQ  between  Thin  Plgte  an^  §SS!S  Theory 
Result^ 


In  chapter  3,  the  critical  buckling  stress  of  a 
flat  bar  stiffener  was  derived  through  the  use  of  thin 
plate  theory: 


2  2  2 

D  d<m>^/a>  Cd  <m77-/a)  +  6<l-V)3/3  +  C  B 

w  o 

t  - -  - 

CR  2  * 

(M  7^/a)  1  B  -t-  B 

p  o  o 


(3-3) 

In  chapter  4,  the  critical  buckling  stress  of  a 
flat  bar  stiffener  was  found  using  one  dimensional  beam 
theory: 


ft 

CR 


2  4  2 

E(1  s  C  )(m7^/a)  -i-  SJCm/T/a)  C  B 


(m*?^/#)  I 


B  -I-  B 
o  o 


(4-3) 
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It  should  be  noted  here  that  when: 


2  2  2 

D  d<m;^/a)  Cd  (ra  TT'/a)  +  6(l-V)3/3  » 

w 

2  4  2 

Ed  s  +  C  XmTr/a)  +  GJ(m77'/a) 
z  w 

the  agreement  between  the  two  methods  would  be  exact. 

Por  the  elastic  range  of  experimentation  and  with: 

3 

Et 

w 

D  » -  G  «  E/2(l+V>  «  E/2.6 

12(1-V) 


I  «  <d  t  )/12 
z  w  w 


(t  •  d 


J  ■  (d  t  )/3 
w  w 

V  -  0.3 


w 


The  above  relationship  reduces  to: 


2  2 
d  (m77'/a) 

w 


2  2 

d  (mr-/a)  1 

w 


10.  92  2. 6 

(Thin  Plate  Theory) 


12  7.8 

(Beam  Theory) 


From  a  cursory  examination  it  would  appear  that  the 
Thin  Plate  Theory  results  (equation  3-3)  will  probably 
give  a  larger  value  for  the  buckling  stress  for  the 
normal  range  of  stiffener  dimensions  used  in  ships.  To 
test  this  hypothesis,  the  dimensions  of  a  flat  bar 


S8 


stiffener  used  in  experiment  #21  of  reference  (6)  will 


be  used, 
where: 

d  »  87.  4  cm 
w 

a  s*  730  cm 
m  ■  1 

For  the  thin  plate  theory: 

2  2 

d  (mTT/a)  1 

_  +  -  ■  0. 398 

10. 92  S<  6 

For  the  beam  theory: 

2  2 

d  <m  TT/a)  1 

w 

_  + - -0. 14 

12  7.8 

Thus,  within  the  elastic  range,  the  thin  plate 
theory  results  would  indeed  provide  higher  buckling 
stress.  The  beam  thoery  results  may  thus  be  too 
conservative.  However,  section  5.3  shows  good  agreement 
between  equation  (4-11)  (first  yield  load)  and 
experimental  yield  results. 


5.2  Cowoariaons  wifeh  Pu^listlSd  Formula^. 


Rafarence  (3)  performed  a  similar  analysis  to  that 
performed  in  this  thesis  but  without  the  inclusion  of 
initial  geometric  imperfections. 

For  the  case  of  a  Tee  stiffener  subjected  to  an 
axial  loadf  equation  21  of  reference  (3)  allows  (in  the 
notation  of  this  thesis): 

2  2  2 
GJ  +  Ed  s  +  C  )(m7^/a)  +  Cia/ta'ii') 

z  w 

CR 

X 

P 


Equation  (2-2)  of  this  thesis,  for  a  Tee  stiffener 
in  an  analysis  which  ignores  the  contribution  of  the 
web  and  includes  the  contributions  of  initial  defects, 
is: 


2  2 

CE(I  s  *  C  6J3  B 


z  w  o 


I  B  >  B 


p  o  o 

Thus  equation  (2-2),  for  the  case  of  an  initially 
perfect  Tee  stiffener  reduces  to  equation  21  of 
reference  (3)  when  the  rotational  restraint  is  held  to 
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p 


be  zero.  However,  it  must  be  noted  that  since  the 

physical  models  of  the  two  derivations  were  slightly 

different  -  namely  that  the  model  used  in  this  analysis 

consisted  of  two  unequal  flanges  some  of  the 

parameters  of  the  equations  will  have  different  values 

(ie.  1  ,  J  etc. ) . 
z 

For  the  case  of  a  flat  bar  stiffener,  equation  22 

of  reference  <3)  has: 

2  2  2 
D  dC(m  ^/a)  d  6(1- V  )3 

w 


C(a/m  /T) 


«  - 

CR 
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From  chapter  4,  equation  (4-3) : 

2  4  2 

Ed  s  C  )  <m /T-Za)  +  BJ(m  r^/a)  +  C 

z  w 


B 


CR 


I  <m  TT'/a) 
P 


« 

B  *  B 
o  o 


The  following  is  a  comparison  for  a  Tee  stiffener 
and  a  flat  bar  stiffener  which  was  presented  in 
reference  <3).  Figure  (5-1)  from  reference  (3)  shows 
the  geometric  dimensions  of  the  Tee  and  flat  bar 
stiffeners.  Table  (5-1)  shows  how  equations  21  and  22 
of  reference  (3)  and  equations  (2-2)  and  (4-3)  of  this 
thesis  compare  with  a  finite  element  analysis  performed 
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1 


using  a  computsr  program  developed  and  documented  at 
the  University  of  California  (reference  12). 


4 

Flat  Bar  (&  /E)  x  10  C  »  0 

CR 

f 

FEM 

E(3N(22) 

EQN(4-3) 

m  a  1 

6.  44 

6.42 

6.  28 

ffl  w  2 

7.sa 

7.51 

7.03 

Tee 

4 

(«  /E)  M  10  C  »  0 

CR 

FEM 

EQN(21) 

EQN(2>2) 

B 

26.4 

26.3 

30.69 

ffl  <■  2 

65.9 

72.6 

89.4 

Table  (S-1)  Comparison  with  Published  Formulae 
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used  similar 


The  U. S.  Navy  has  in  the  past 
formula  as  those  presented  here  for  the  determination 
of  critical  buckling  stresses.  For  the  critical 
buckling  stress  of  a  Tee  stiffener,  in  the  notation  of 
this  thesis  from  reference  13  Me  have: 

2 

6J  *  (mTT/a)  EC 

M 

ft 

CR 

I 

P 

where  here: 

3  2  3  3 

C  «  f  d  t  /12  +  d  t 
.w  w  f  w 

3  2  3 

I  -dt/3  +  dft  ■►ft 
p  w  w  w  w  f 


This  is  very  similar  to  those  results  determined  by 
reference  (3),  with  the  exception  that  the  above 
formula  does  not  account  for  the  rotational  restraint 
provided  by  the  junction  of  the  web  and  plate.  Also  I 

P 

is  much  simplified  here  and  C  ,  the  torsion-bending 

w  2 

constant,  equates  to  the  quantity  (Is  -►  C  )  of 

z 

reference  (3)  -  but  also  somewhat  simplified. 


6A 


w 


5.3  Comoariaons  wi.th  E x pgr i mental  Results. 

Reference  (&)  provides  insight  into  the  torsional 
buckling  process  through  the  experimental  results 
published  there.  This  report  describes  model  tests  on 
the  collapse  and  post  failure  strength  of  Tee  struts 
simulating  flat  bar  stiffeners  in  a  stiffened  plate. 
The  T  struts  were  subjected  to  an  axial  load  up  to  and 
beyond  failure.  The  geometrical  imperfections  of  the 
models  were  also  recorded.  The  researchers  of  reference 
(6)  (of  Det  norske  Veritas)  are  to  be  lauded  for  the 
thoroughness  and  completeness  of  their  report, 
particularly  in  the  treatment  of  the  experimental  data. 

Table  <5-2)  contains  the  model  data  pertaining  to 
this  analysis  for  ten  experiments  of  reference  (8). 
Figures  (5-2  to  5-11)  show  the  load-  deflection  curves 
for  these  ten  experiments.  The  solid  line  represents 
the  experimental  data.  Point  A  of  each  of  these  graphs 
represents  the  first  yield  load  point  as  determined  by 
equation  (4-11). 

Notice  that  there  is  very  good  correlation  between 
the  experimental  first  yield  load  point  and  the  results 
of  equation  (4-11). 
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Model  a 

h 

t 

B 

ft 

ft 

ft 

w 

o 

y 

CR 

u 

21 

703 

87.4 

4.0 

1.67E-02 

256 

233 

230 

22 

703 

87.3 

4.0 

5. 6E-03 

271 

263 

246 

23 

703 

87.7 

4.0 

1. lSE-02 

273 

257 

234 

24 

703 

87.7 

4.0 

6. 0E-03 

271 

262 

240 

25 

703 

87.7 

4.0 

8. 0E-03 

272 

260 

230 

26 

703 

87.7 

4.0 

8. 0E-03 

256 

244 

243 

27 

703 

87.4 

4.0 

6. 0E>03 

266 

257 

272 

28 

703 

88.0 

4.0 

7. 99E-03 

256 

244 

241 

29 

703 

87.8 

4.0 

4. 96E-03 

259 

252 

218 

30 

703 

88.  1 

4.0 

1.26E-02 

271 

253 

230 

Notation: 

I 

a  ■*  length  of  model  (cm) 

h  -  height  of  stiffener  (cm) 

t  -  thickness  of  stiffener  (cm) 
w 

* 

B  -  angle  of  rotation  (radians) 
o 

2 

ft  -  yield  stress  (N/mm  ) 

y 

2 

ft  -  first  yield  load  stress  (N/mm  ) 

CR 

2 

ft  -  ultimate  stress  (sMperimental  results) (N/mro  ) 
u 


Table  (S-B)  Data  of  Model  Tests 
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STR0<ECW1>  t  M* 

MODEL  2t 


Figur*  (5-3) 


Mod»l  TMt  #2^  of  Roforonco  (6) 


•  o.oct  «.70 
MODEL  2$ 


1.40  2.10  2.80  3.50  •4.C0  4.90  S.bg  5.3« 

stroke «nn> 


Figur*  <5-6)  Modal  T»»t  #85  of  Rafaranca  (6> 
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EXPCRIMENt 


HOOEl  26 


Fiflur*  (5-7)  Mod«l  T«»t  »S6  of  Roforonco  (6) 


72 

'■  n 

1 


«.7e  1.49  8.10  8.89  3.59  4.89  4.50  S.eO  6. 3.»  7 


STSOKClIirt)  L 

MODEL  28 


Flgur«  (5-9)  VMml  T«»t  #20  of  R«f»r«nc«  (6) 


STR0IC£<|1I1>  C  !»■» 


MODEL  29 


Figur*  (S*-19)  Modal  Tast  #29  of  Rafaranca  (6> 
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Figure  (&-11)  Model  Teet  #30  of  Reference  (6> 


6.  CONCLUSIONS 


This  thesis  utilized  a  design  approach  to  tackle 
the  problem  of  torsional  buckling  of  stiffeners  with  a 
good  deal  of  success.  Integrated  into  all  phases  of 
this  analysis  was  a  consideration  of  the  geometric 
imperfections  of  the  stiffener  -  both  Tee  and  flat  bar. 
These  structural  members  were  analyzed  as  simplified 
sections  and  the  Perry-Robertson  approach  used  in  other 
areas  of  engineering  was  successfully  applied  here. 

Chapter  S  showed  that  the  derived  results  of  this 
thesis  reduced  to  currently  published  formulae  for  the 
case  of  the  perfect  stiffener.  The  contribution  of  this 
thesis  is  that  the  initial  imperfections  were  taken 
into  consideration,  using  a  simple  design  approach,  in 
the  derivation  of  the  results  presented  here. 

Also  shown  in  chapter  5  was  a  generally  good 
correlation  between  results  of  this  thesis  and  axial 
load  tests  on  flat  bar  stiffeners.  Comparison  of  the 
Tee  stiffener  formula  with  experimental  data,  which 
measures  and  reports  initial  imperfections,  needs  to  be 
performed.  Also,  comparisons  should  be  made  with  other 
tests  on  flat  bar  stiffeners  in  more  dimensional  ranges 
to  determine  the  range  of  validity  for  the  results 
derived  here. 


Further  work  needs  to  be  done  in  this  area. 


Specifically,  determination  of  the  critical  buckling 
stress  and  first  yield  load  due  to  lateral  loadings. 
Also,  a  simplified  method  of  handling  the  mode 
interactions  of  combined  lateral  and  axial  loadings  is 
needed.  In  addition,  it  is  an  asset  to  understand  from 
the  work  reported  in  references  (S)  and  (15)  that  the 
falling  path  of  the  load-def lection  curve  (figure  &-1) 
is  represented  by  the  tripping  mechanism  which  is  not 
covered  in  this  analysis. 


!•  I-  a  U  'A 


E 


MOOEl  n 


Figure  (6*^1)  EMperiiMintel  Load  Deflection  Curve 
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